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Abstract 

A novel unitary transformation of the Hamiltonian that allows one to par- 
tially separate the center-of-mass motion for charged electron-hole systems in 
a magnetic field is presented. The two- mode squeezed oscillator states that 
appear at the intermediate stage of the transformation are used for construct- 
ing a trial wave function of a two-dimensional (2D) charged magnetoexciton. 
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A problem of center-of-mass (CM) separation for a quantum-mechanical system of 
charged interacting particles in a magnetic field B has been studied by many authors.^! 
When a charge-to-mass ratio is the same for all particles, the CM and internal motions 
decoupleEHl in B. For a neutral system, the CM coordinates can be separated^! in the 
Schrodinger equation. This is associated with the fact that translations commute for a 
neutral system in B. In general, only a partial separation of the CM in magnetic fields is 
possibleBBi In this Letter we propose a novel operator approach for performing such a sepa- 
ration in charged electron- hole (e-h) systems in B. This approach can be useful for studying 
in strong magnetic fields, e.g., atomic ions with not too large mass ratios! and charged ex- 
citations in 2D electron systems, in particular, in the fractional quantum Hall effect regime 
in the planar geometry.!0 In this work, we study in 2D a three-particle problem of two 
electrons and one hole in a strong magnetic field, i.e., a negatively charged magnetoexciton 
X~ (see Refs. HHH and references therein). We consider an approximate X~ ground state 
in the form, which is related to the two-mode squeezed^ oscillator vacuum states. 

The Hamiltonian describing the 2D three-particle 2e— h complex in a perpendicular mag- 
netic field B is H = Ho + H mt , where the free-particle part is given by 

TT 2 TT 2 

^ 2m e 2m h ^ 0eK l) 0hK h> ' K J 

and Uj = —ih'Vj — — A(r^) are kinematic momentum operators. The interaction Hamilto- 
nian H int = H ee + H eh is 

e 2 e 2 

H ee = -j r 7 H ch = — 2J -j r • (2) 

e| r i — r 2| j = i 2 e| r i — r h\ 



The Hamiltonian H commutes™™ with the operator of magnetic translations (MT) K = 
£),-Kj, where Kj = IT,- — x B. In the symmetric gauge, A = |B x r, the operators 
satisfy the relation K,(B) = ILj(— B); independent of the gauge, Kj and flj commute. The 
important feature of K and n = J2j LL, is the non-commutativity of the components in B: 
[K x ,Ky] = —[fl x ,fly] = —i^f-Q, where Q = J2j e j is the total charge. This allows one to 
introduce the raising and lowering Bose ladder operators for the whole systemoQ'B 

k± = ±—^=(k x ± iky) , [k + , k-] = — , (3) 
i Q 



where k = ^Jc/hB\Q\K., fr = Jc/hB\Q\tl, and the phases of the operators (|3|) and (f|) 

can be chosen arbitrary. The operator k 2 has the discrete oscillator eigenvalues 2k + 1, 
k = 0, 1, . . . that are associated!!! with the guiding center of a charged complex in B. The 
values of k can be used,0 together with the total angular momentum projection M z and the 
electron, S e , and hole, Sh, spin quantum numbers, for the classification of states;! the exact 
eigenenergies are degenerate!!! in k. 

In terms of the single-particle Bose ladder intra Landau level (LL) operators^! B\{rj) = 

—iJc/2hBe (Kj x — iKjy) f° r the electrons and B^r^) = —iJc/2hBe (Khx + iKhy) f° r the 
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hole, the raising operator takes the form k- = B\{ji) + Bl(r 2 ) — Bhijh). One needs to 
diagonalize /c_ in order to maintain the exact MT symmetry. This can be achieved by 
performing first an orthogonal transformationll! of the electron coordinates {r 1 ,r 2 ,r/ l } — > 
{r, R, Vh}, where r = (ri — r 2 ) / y/2, and R = (ri + r 2 ) / y/2 are the electron relative and CM 
coordinates. In these coordinates fc_ = y/2B\{R) — Bh(rh) and can be considered to be a 
new Bose ladder operator generated by the Bogoliubov transformation! 



Bt(R) = uB\{K) - vB h (r h ) = SBt(R)& 



(5) 

where the unitary operatorEHiiJ S = exp(0£) and the generator £ = S^(r^)5|(R) — 
B e (R)Bh(rh). Here O is the transformation angle and u = cosh0 = y/2, v = sinh0 = 1. 
Now we have = B\ and k = 2B\B e + 1. The second linearly independent creation 
operator is 



Bl(v h ) = SBl(r h )& = uBl(v h ) 



vB e (R) 



(6) 



A complete orthonormal basis compatible with both axial and translational symmetries can 
be constructed! as: 



At(r)"^t(R)"^j(r fe )"^t(R)^t(r fe )'i?t(r) m |0) 
(rirlnjilnhlkHlml) 1 ^ 2 
= \n r n R n h ; klm) . 



(7) 



Here the inter-LL Bose ladder operators are given by A\{vj) = — c/2hBe (tlj X + itlj y ) and 
Al(rh) = — i\J~c/2hBe (flh x ~ i^-hy)] the explicit form is given in, e.g., Refs. |||]. The tilde 
sign shows that the transformed vacuum state |0) (see below) and the transformed operators 
(I) and @ are involved. In (0) the oscillator quantum number is fixed and equals k, while 
M z = n r + tir — rih — k + I — m. The permutational symmetry requires that n r — m should 
be even (odd) for electron spin-singlet = (triplet S e = 1 states); see Ref. || for more 
details. 

The transformation introduces a new vacuum state |0) = S^O), for which, using the 
normal-ordered form!Ji0 of S, one obtains 



|0> = S\0) 



exp 



cosh G 

The coordinate representation has the form 



tanh 6 B\(r h )Bl("R) |0) . 



(rRr ft |6) 



1 



V2(2nl 2 B ) 3 / 2 



exp 



r 2 + R 2 + r 2 



y/2Z* 



All 



(8) 



(9) 



where Ib = (hc/eB) 1 ^ 2 is the magnetic length, Z* = X — %Y , and Zh = Xh + iyh- Equation 
(§) shows that |0) contains a coherent superposition of an infinite number of e- and h- states 
in zero LL's. In the terminology of quantum optics," |0) is a two-mode squeezed state; for 
particles in a magnetic field the squeezing has a direct geometrical meaning.0 Indeed, the 
probability distribution function takes the factored form 
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8/| 



;R+r„) 2 



(10) 



This shows that the distribution for the relative coordinate R — is squeezed at the expense 
of that for the coordinate R + r^, and the variances are (0|(R ± r^) 2 10) = 4(2 ± \/2) l\. 
The squeezing enhances the e—h attraction which will be used below for constructing a trial 
wave function of the 2D magnetoexciton X~ . 

Let us now perform the second unitary transformation corresponding to the diago- 
nalization of the operator^ 7r + = A\(ri) + A\(r 2 ) — Ah(rh). This introduces a new 
state |0) = 55 1 0) = 5|0), which corresponds to the simultaneous diagonalization of 
the operators k_ and 7r + ; the unitary operator 5 = exp(0£), where the generator 
L = A^ThjAlCR) — A e (R)A h (r h ). The transformations effectively introduce new coor- 
dinates, {r, R, r^} — > {r, pi, p 2 }, where pi = \/2R — and p 2 = V^Th — R, which can be 
presented in the matrix form 





cosh O — sinh B 
— sinh O cosh B 



(11) 



with coshO = y/2, sinh0 = 1; the matrix F corresponds to the SU (1,1) symmetry.0 
Indeed, the inter-LL ladder operators are changed under the Bogoliubov transformations as 

S {A h (r h )) S - F \A h (v h ))-\A h (p 2 )) ■ W 

The intra-LL operators ([5]) and (§) transform according to the same representation. The 
coordinate representation 

<W 2 |0H^exp(-£!±||±ii) (13) 

shows that |0) is a true vacuum for both the intra-LL B\(pi), Bj l (p 2 ) and inter-LL A\(p 2 ), 
A\(pi) operators. Now we can perform the change of the variables {r, R, r^} — > {r, pi, p 2 } 
in the basis states: 

\n T nnnh'i klm) = 

(nrlnnlrihlkHlml) 1 / 2 
= S^\n r n R n h ; klm) . (14) 

The overline shows that a state is generated in the usual way by the intra- and inter-LL Bose 
ladder operators acting on the true vacuum |0) — all in the representation of the coordinates 
{r, pi,p 2 }. The Hamiltonian H is block-diagonal in the quantum numbers k,M z (and 5 e , 
Sh)- Due to the Landau degeneracyi'i in k, it is sufficient to consider the states with k = 0. 
This effectively removes one degree of freedom and corresponds to a partial separation of the 
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CM motion. From now on we will consider the k = states only, designating such states in 
(14) as \n r n R nh] lm). For the Hamiltonian we arrive therefore at the unitary transformation 

(m 2 l 2 ; n h2 n R 2n r2 \H\n rl n R1 n hl ; l x m\) 

= (7712/2; n h2 n R 2n r2 \SHS ] \n rl n R1 n h i] l x mi) , (15) 

which is the main formal result of this work. 

The Coulomb interparticle interactions (0) in the coordinates {r, pi, p 2 } take the form 

H ee = f- , H eh = —— 7^— , (16) 

V2er e\p 2 - r\ e|p 2 + r| 

and H mt does not depend on p x . From Eq. fll2|) it follows that the free Hamiltonians trans- 
form as SH 0e (r)S^ = H 0e (r), SH 0e (R)S^ = H Qe (pi), and SH oh (r h )S^ = H 0h (p 2 ) and de- 
scribe new effective particles — free e and h in a magnetic field — with the modified inter- 
actions fllBl) . The Hamiltonian of the e-e interactions H ee (\/2\r\) does not depend on pi, p 2 
and, therefore, is invariant: SH ee S^ = H ee . Thus, the matrix elements of the e-e interaction 
are easily obtained from ([Tj]): they reduce to the matrix elements V n ymx describing the 
interaction of the electron with a fixed negative charge — e: 

(m 2 l 2 ; n h2 n R2 n r2 \H ee \n r in m n h i; l x mi) (17) 
= (m 2 / 2 ; n h2 n R2 n r2 \H ee \n r in R in h i; hmi) 

_ X X A A 1 T/n r 2 ni2 

u n R1 ,n R 2 u n hl ,n h 2 U ll,h u n r i-m 1 ,n r 2-m 2 ^ 2 v n rl m 1 • 



In, e.g., zeroQ LL Vq™ = [(2m — l)\\/2 m m\] E , where E = «/| J^-. The generator C and 
the Hamiltonian H e h(r, p 2 ) do not form a closed algebra of a finite order. Therefore, the 
explicit form of SH^S* cannot be found. We can find, however, the form of the matrix 
elements of SH^S* in (|T5|). Because of the electron permutational symmetry r <->• — r it is 
sufficient to consider the term U e h(p2 — r) = —e 2 /e\p 2 — r|. Here we only consider the states 
in zero LL 1 000; lm) = \lm). Using the normal-ordered form of S, we have 

(7nA\SU eh &\hm^) ee (18) 
= |(^^|e-^^ {pi) ^ (p2 V eh e-^ Al(p2Mi(pi) |77^) . 

Expanding the exponents and exploiting the fact that U c \ l (p 2 — r) does not depend on pi, 
we obtain a series 

00 

fjQrm Oh _ 1 ( l\ p TT0m 2 pl2 (-1 q\ 

^OmiOh ~ 2 Z^\2J U 0m lP h ■ \ LJ ) 
p=Q 

Note that ([19]) includes contributions of the infinitely many LL's. For the Coulomb inter- 
actions the matrix elements can be calculated analytically; in zero LL we obtain 

(m 2 l 2 \SH eh S Jf \limi) = 8 h 

—m\,l2-m2 

2\/2 ?7niin(mi,m2),min(Ji,J2)(l m l ~ 77l 2 |) , (20) 
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[m\(m + s)\n\(n + s)!] 1/2 2 m + n + s 3 s+1 / 2 
xEE^X (!) m + 1 + 8)- 1]!! [2(m - k) - 1]!! 

fc=0 Z=0 

xE^Cf-lf?! [2(n-J-p)-l]H. 

p=0 

The developed formalism can be used for performing a rapidly converg enti expansion of 
the interacting e—h states in the basis dl5|), which preserves all symmetries of the problem. 
Here we demonstrate a possibility of using the squeezed states for constructing trial wave 
functions. We consider the triplet charged 2D magnetoexciton in zero LL, X^ , with M z = 
— 1, which is the only bound state!'! in zero LL in the strictly-2D system in the high-field 
limit. The simplest possible wave function in zero LL compatible with all symmetries of the 
problem is 

This form allows analytic calculations and, as a squeezed state (see above), already 
ensures the X 4 g binding. Indeed, the total Coulomb interaction energy is given by 
-^oV + 2^^10(0) = - ^-)E ~ -l.007.Eo. The corresponding binding energy 
(counted from the ground state energy of the neutral magnetoexciton, —E ) is 0.007-E , 
which is 17% of the numerically "exact" value of 0.043i?oB! A similar type of squeezing 
can be applied to construct a trial wave function of the X^ . The idea is to additionally 
squeeze the effective hole p 2 and electron r coordinates. Since the wave function must 
by antisymmetric under the permutation of the electron coordinates, we can use the form 
|0) ~ Bl(r)(S ( j > + S-^)S^\0), where the second two-mode squeezing operator is given by 
5*0 = exp[(j)Bl(r)Bj l (p2) — H.c] and we have used |0) = «S^|o). The normalized four-mode 
squeezed wave function has the form 

l + tiUllrn £t( r ) cosh\tanh<f>Bl(r)Bl{p 2 )]S^0). (23) 



cosh 2 <b\ll + tanh 



The calculated energy of the Coulomb e-e repulsion, (0|iJ ce |0), monotonically decreases with 
increasing the transformation angle 0, whereas the energy of the e—h attraction, — (0|iJ e h|0), 
has a maximum (see Fig. 1). The binding of the X 4 q results from a rather delicate balance 
between the two terms, and for the state (|23| ) the maximum achieved binding energy is 
E b ~ 0.038i?o (see inset to Fig. 1), which is 91% of the "exact" value;!'! note that the 
inaccuracy is 0.3% of the e—h interaction energy. Similar type of squeezed trial wave functions 
may be useful in other solid state and atomic physics problems dealing with correlated e—h 
states in strong magnetic fields. 

In conclusion, we have developed for charged e—h systems in magnetic fields an operator 
approach that allows one to partially separate the CM motion. This results in the appearance 
of new effective particles, electrons and holes in a magnetic field, with modified interparticle 
interactions. A relation of the considered basis states with the two-mode squeezed oscillator 
states has been established. 
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FIGURES 

FIG. 1. The expectation values of the e-e repulsion, ((p\H cc \<p) , and the e—h attraction, 
(4>\H c h\(f)) (with the opposite sign, counted from the neutral magnetoexciton binding energy 
Eq = f° r ^ ne trial wave function ( |23|) of the charged triplet magnetoexciton in zero LL's, 

X t Q . The binding energy E\, = — (<^|H e h + ^ee|<^) ~~ ^0 is shown in the inset. 
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Transformation angle § 



